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ABSTRACT 

The W28 supernova remnant is an excellent prototype for observing shocked gas 
resulting from the interaction of supernova remnants (SNRs) and adjacent molecular 
clouds (MCs). We present two new signatures of shocked molecular gas in this 
remnant. One is the detection of main-line extended OH (1667 MHz) absorption with 
broad linewidths. The column density of OH estimated from the optical depth profiles 
is consistent with a theoretical model in which OH is formed behind a C-type shock 
front. The second is the detection of extended, weak OH (1720 MHz) line emission 
with narrow linewidth distributed throughout the shocked region of W28. These give 
observational support to the idea that compact maser sources delineate the brightest 
component of a much larger region of main line OH absorption and nonthermal OH 
(1720 MHz) emission tracing the global structure of shocked molecular gas. Main line 
OH (1665/67) absorption and extended OH (1720 MHz) emission line studies can serve 
as powerful tools to detect SNR-MC interaction even when bright OH (1720 MHz) 
masers are absent. 

Subject headings: ISM: Clouds — ISM: general — shock waves — supernova remnants — X- 
rays: ISM 
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1. Introduction 

OH (1720 MHz) masers unaccompanied by masing in the 1665 and 1667 MHz main lines 
occur where supernova remnants (SNRs) interact with molecular clouds (MCs) (Frail, Goss &; 
Slysh 1994; Yusef-Zadeh et al. 1996; Green et al. 1997; Wardle and Yusef-Zadeh 2002). These 
so-called SNR masers arc coUisionally pumped within a shock wave driven into a dense molecular 
cloud. The OH molecule is believed to be produced by the indirect dissociation of water formed 
within a C-type shock (Lockett, Gauthier k, Elitzur 1999) by X-ray photons from the interior of 
the supernova remnant (Wardle 1999). 

Despite several radio surveys to find SNR-MC interaction sites in the Galaxy, OH (1720 MHz) 
masers have been detected toward only 19 SNRs, about 10% of the total number of observed SNRs 
(Frail et al. 1997; Koralesky et al. 1998; Green et al. 1997; Yusef-Zadeh et al. 2002). This low 
detection rate is partly due to the fact that OH (1720 MHz) masers only arise under restrictive 
conditions in density and temperature (Lockett et al. 1999) and that velocity coherence and 
maximum pathlength for amplification are achieved at the limb of the remnant where the shock 
propagates in a direction perpendicular to the line of sight. However, the OH abundance increases 
throughout the shocked region and is expected to be more extensive and less restrictive in its 
physical conditions (Wardle 1999) than the regions where bright maser spots occur. 

The OH column densities required for the formation of OH (1720 MHz) masers 10^^ cm~^) 
should be sufficient to create detectable thermal OH absorption in the main line transitions of OH 
(1665/67 MHz) molecule if the background continuum is strong. In order to test this idea, we 
observed a nearby interacting SNR, W28 at the distance of 1.8 kpc (Rho et al. 1994). W28 is a 
prototype SNR maser with numerous compact and bright maser spots distributed throughout the 
remnant (Claussen et al. 1997). OH absorption towards W28 has previously been detected with 
single-dish observations (Pastchenko & Slysh 1974) and with the Australian Telescope Compact 
Array (Green et al. 2000). Here we report the VLA detection of broad OH (1667 MHz) absorption 
lines against the prominent W28 SNR continuum source. The column density of shocked OH gas 
estimated from optical depth profiles toward W28 is consistent with the theoretical prediction. In 
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addition, we present extended, weak OH (1720 MHz) maser emission from W28. We then suggest 
that studies of thermal OH (1667 MHz) absorption and extended nonthermal OH (1720 MHz) 
emission can be used to search for additional SNR-MC interactions where compact 1720 MHz 
maser s are absent. 



2. Observations 

The Very Large Array (VLA) of the National Radio Astronomy Observatory^ was used on 
October 7, 2001 and July 18, 2000 to observe OH (1667 MHz) absorption and OH(1720 MHz) 
emission line studies, respectively. We used the compact D configuration to observe OH (1667 
MHz) line absorption toward W28 at center velocities of 10 kms~^. Mode 4 of the correlator was 
used with 64 channels per IF; each IF pair is tuned to the different line transitions of OH, (1667 
MHz and 1720 MHz). The channel separation was 6.104 kHz, which corresponds to 0.55 km/s. 
Standard calibrations of complex gains and bandpass were carried out using 1328+307, 1730-130 
and 1851-235. Because of a strong negative bowl produced surrounding the continuum image of 
W28, a zero-spacing flux of 60 Jy was added to the continuum uv data. The uv data were then 
self-calibrated before the final line (1) and continuum (c) images were CLEANed. The rms noise in 
a given channel of W28 is 5.5 mJy beam~^. An optical depth (r) image and corresponding error 
images were constructed for W28 by using r = —ln{{l + c)/c). 

In the second VLA observation, we carried out an OH (1720 MHz) emission line study of 
W28. Unlike the previous high-resolution snapshot OH (1720 MHz) observations of W28 using 
the A-configuration of the VLA (Frail et al. 1994; Claussen et al 1997), we searched for extended, 
weak OH (1720 MHz) emission from W28 using the compact D configuration with a relatively 
uniform uv coverage so as to be sensitive to extended but narrow line OH (1720 MHz) emission 
features. Mode 2AC of the correlator was chosen with 128 channels, a channel width of 1.06 kms~^ 

^The National Radio Astronomy Observatory is a facility of the National Science Foundation, 
operated under a cooperative agreement by Associated Universities, Inc. 
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and a total bandwidth of 0.7812 MHz. This observation used similar phase and flux calibrators to 
that of absorption line observations. The rms noise for a given channel is 6.2 mJy beam"^. The 
rms noise for three velocity channels at 9.6, 10.7 and 11.8 kms~^ where the most intense maser 
emission is observed were three times higher. 

3. Results 

3.1. W28 (SNR G6.4-0.1) 

3.1.1. Broad Thermal OH (1667 MHz) Absorption 

The right panel of Figure 1 shows a contour representation of OH (1667 MHz) line absorption 
toward W28 integrated between -3.2 and 23.7 kms~^. The grayscale image shows the continuum 
image of W28 constructed from line- free channels. Much of the absorption is distributed 
prominently to the northeast of the remnant where background continuum is strongest. The 
absorbing gas velocities are very similar to OH (1720 MHz) maser velocities ranging between +4.8 
and +16 kms~^ as reported by Claussen et al. (1997) and Frail et al. (1994). The panel to the 
left of the remnant in Figure 1 shows a typical absorption spectrum (position is indicated by the 
cross on the right panel) with a total linewidth of «30 kms~^. 

Figure 2 shows a grayscale continuum image at 20cm (Dubner et al. 2000) with the 
distribution of compact OH (1720 MHz) maser spots represented as crosses (Claussen et al. 1997) 
superposed. The insets to the left and top of Figure 2 show the fitted optical depth profiles and 
their corresponding residuals. The optical depth velocity profiles coincide with clusters of OH 
(1720 MHz) maser sources identified by Claussen et al. (1997). Gaussian fitted spectra of the top 
and left r profiles show center velocities between \lsr = 10-2 and 8.1 kms~^, which are similar to 
the velocities of of OH (1720 MHz) maser spots identified in the F, C and D regions of Claussen et 
al. (1997). The linewidths (FWHM) of 14 and 11 kms~^ are much larger than typical linewidths 
of ambient molecular clouds, thus confirming the interaction picture of W28 with its surrounding 
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molecular cloud (e.g. Wootten 1981). Table 1 gives the parameters of the fits to four r spectra, 
all of which show broad line widths ranging between 7 and 10 kms~^. However, the top-right 
spectrum of Figure 2 shows two narrow absorption lines with r=0.23 and 0.28 plus a broad 
absorption feature with r=0.08 at centroid velocity of 10 kms~^. This spectrum is seen against 
a bright compact continuum source G6.6-0.1 a, 5 (J2000)= 18^00™49.8^ -23°20'24".56 showing a 
flux density of 0.76 Jy near the center of W28. This source is known to have a flat spectral index 
(Andrews et al. 1983; Dubner et al. 2000). The line profiles show FWHM linewidths of 1.5-1.8 
kms^^ and are centered at 12 and 18 kms^^. The narrow line widths of absorption features 
suggest that they are not associated with the SNR-MC interaction. The high velocity component 
is close to the high velocity component of HI absorption measurements (Radhakrishnan et al. 
1972). Claussen et al. (1997) suggest that the systemic velocity of the molecular cloud associated 
with W28 is at 17.6 kms~^ and that the low-velocity 7 kms~^ component is being accelerated by 
the SNR shock along the line of sight by about -10 kms^^. This suggests that the narrow line 
spectra observed in Figure 2 are associated with the MC in the vicinity of W28. 

A bright, compact OH (1667 MHz) maser source with two velocity components is detected in 
the vicinity of W28. This source lies at a, 5 (J2000)= 18''01™57.78", -23°12'31".43. The spectrum 
shows two narrow velocity components at 17.1 kms~^ and 23.5 kms~^ with a flux density of 0.39 
and 0.05 Jy, respectively, resembling the spectrum of OH/IR stars. 

3.1.2. Extended Maser Emission 

Figure 3 shows contours of the OH(1720 MHz) line emission integrated between 5.3 and 21.3 
km/s (with the exception of three channels at 9.6, 10.6 and 11.7 kms~^, which were contaminated 
by saturation from strong maser emission) superposed on a grayscale continuum image of W28. 
The crosses indicate the positions of 41 compact maser sources listed by Claussen et al. (1997). 
We note extended and compact OH (1720 MHz) emission from W28. The top two insets show 
the velocity profiles of two isolated sources to the north of the shell. These sources, as well as 
extended maser emission adjacent to compact maser sources, have not been detected in earlier 
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high resolution observations which over-resolved these extended features. The third OH (1720 
MHz) velocity profile in the bottom left panel of Figure 3 coincides with the position of the optical 
depth profile presented in the top left panel of Figure 2. It is clear that where OH(1667 MHz) 
absorption profile is broad, OH (1720 MHz) line emission is narrow. 

The extended OH (1720 MHz) emission features have typical linewidths of about 1-3 kms~^. 
The flux density of extended features arc typically about 100-300 mJy/beam corresponding to 
a brightness temperature T^ ~ 20-60 K. The kinematics of extended OH (1720 MHz) emission 
is very similar to that of compact masers, suggesting that they co-exist with each other. The 
extended features with their nonthermal characteristics trace shocked CO (3-2) molecular gas 
toward W28 (Arikawa et al. 1999) in spite of their low brightness temperatures. 



4. Discussion 

The column density of OH toward W28 may be estimated from the observed optical depth 
velocity profiles using 

iVoH = 2.2785 X 10^"^ Tex j dv cm'^ (1) 

for the 1667 MHz line once a value for the excitation temperature, T^^, has been chosen (Crutcher 
1977). Because the energy spacing between the upper and lower levels is RiO.08 K, Tex is sensitive 
to small departures of the level populations from LTE. In order to estimate Tex, we modelled OH 
level populations in a uniform slab of gas with a FWHM linewidth of 10 kms~^. Tex is determined 
largely by the kinetic temperature T^ and is insensitive to density and column density for the 
parameters that are relevant here. The excitation temperature is close to the kinetic temperature 
Tfe for Tfc ;^ 20K and then decreases with increasing kinetic temperature as collisions begin to 
populate the higher rotational states: Tex = 18, 8, and 3 K at T^ = 30, 50 K and 75 K respectively. 
Therefore we adopt T^x = 10 K as a nominal value with the assumption that T^x may differ by a 
factor of two in either direction. Our measurements indicate that f dv ^ 10 kms~^, therefore 
Nqu pa 2 X 10^^ cm~^. This OH column density is consistent with that required to produce 1720 



-8- 



MHz maser emission (Lockett, Gauthier & Elitzur 1999) and expected from X-ray dissociation of 
water produced in a C-type shock wave driven into the molecular cloud by the supernova remnant 
(Wardle 1999). 

Extended 1720 MHz emission can be produced by weak amplification in the OH behind a 
face-on shock front as opposed to the edge-on geometry required by the bright maser spots. The 
brightness temperature of the 1720 MHz emission, ~ 30 K at line center, is roughly equal to 
that of the background continuum. The implied amplification factor of 2 requires an OH column 
NoH ~ lO^^cm-2 with ;^ 50K (Lockett et al 1999), and line FWHM ^ 0.5kms-\ conditions 
that are achieved in the cooling tail of a C-type shock wave (Wardle 1999). It should be possible 
to detect this weakly-masing component in absorption at 1667 MHz, as the optical depth is 2 
at line center. H2CO absorption study of W28 also indicates that the interacting cloud has a 
low kinetic temperature and that the properties of the interacting cloud are consistent with the 
picture of cold post-shocked gas (Slysh et al. 1980). 

The observations presented here trace extended regions of shocked molecular gas in which the 
abundance of OH has been enhanced. Our results suggest two new ways to search for SNR-MC 
interaction sites where OH (1720 MHz) masers are absent. One is to search for extended OH (1720 
MHz) maser emission produced by shock waves. Unlike bright, compact masers this emission does 
not require the shock to propagate almost perpendicular to the line of sight, and should therefore 
be more generic. This emission has been detected in the five SNRs with compact masers that 
have been searched to date (Yusef-Zadeh et al. 1995; 1999; and this paper). SNR-molecular cloud 
interaction sites can also be identified by searching for broad OH absorption lines at 1665/67 MHz 
delineating the region where shocked molecular gas overlies the radio continuum emission from the 
SNR. A strong empirical association between X-ray emitting centrally-filled mixed-morphology 
remnants and mascr-cmitting shell-type SNRs has recently been demonstrated by Yusef-Zadeh et 
al. (2002). Remnants that show centrally filled thermal X-ray emission are strong candidates for 
such studies. 

Acknowledgments: We thank G. Dubner for providing us with the continuum image of W28. 
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Fig. 1. — [right panel] Contours of OH (1667 MHz) absorption toward W28 integrated between 
—2.1 and 23.7 kms~^ at —4.0, —2.0, —1.0, —0.5 mJy beam^^ kms~^ are superposed on a grayscale 
continuum image with a spatial resolution of 68" x 30", PA=40°. The range is grayscale brightness 
is between and 270 mJy. [left panel] A typical spectrum showing an OH (1667 MHz) absorption 
line toward the cross at a,5 (J2000)= 18'^01"^41.43% -23°25'11".3. 



Fig. 2. — A grayscale continuum image at 20cm with a spatial resolution of 88" x 48" (Dubner 
et al. 2000) with the distribution of compact OH (1720 MHz) maser spots represented as crosses 
(Claussen et al. 1997) superposed. The insets show the derived optical depths (including error 
bars), Gaussian fits, and corresponding residuals. 



Fig. 3. — A grayscale continuum image of W28 at 20cm (Dubner et al. 2000) with superposed 
contours of OH (1720 MHz) maser emission set at -3, 3, 5, 7, 9, 11, 13, 15, 20, 25, 30, 40, 50, 
60, 70, 80, 100 mJy beam"! with a spatial resolution of 86.6" x 34.6", PA=14° (shown as the 
white ellipse in the bottom right hand corner). The crosses represent the position of compact 
OH (1720 MHz) masers identified by Claussen et al. (1997). The respective positions of the 
three spectra on the top, bottom left, and top left are: a, 5 (J2000)= 18''01°'43.6", -23°24'19".8, 
18^01"'43.33",-23°12'59".2 and 18''01"'32.4^ -23°ll'l0".01. 
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Table 1. OH (1667 MHz) Optical Depth towards W28 



Position 


Amplitude 


Center Velocity 


FWHM Line Width 


Corresponding OH 


(RA, Dec (J2000)) 


(optical depth) 


(km s~^) 


(km s~^) 


(1720 MHz) masers 


18''00"^49?924, -23°20'26^'66 


0.229 ± 0.017 


12.471 ± 0.130 


4.439 ± 0.387 


A 




0.276 ± 0.018 


18.870 ± 0.093 


3.328 ± 0.253 






0.085 ± 0.011 


9.983 ± 1.156 


23.980 ± 2.687 




18''01™50!002, -23°18'46'/58 


0.726 ± 0.027 


10.269 ± 0.183 


14.163 ± 0.765 


F 


18''01°^44!041, -23°24'23?05 


1.148 ± 0.075 


6.845 ± 0.195 


8.928 ± 0.489 


C, D 


18''01°^41!444, -23°25'05?26 


0.668 ± 0.026 


8.533 ± 0.170 


11.508 ± 0.532 


C, D 
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